Interplay of electromagnetic noise and Kondo effect in quantum dots 
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We investigate the influence of an electromagnetic environment, characterized by a finite 
impedance Z(u>), on the Kondo effect in quantum dots. The circuit voltage fluctuations couple 
to charge fluctuations in the dot and influence the spin exchange processes transferring charge be- 
tween the electrodes. We discuss how the low-energy properties of a Kondo quantum dot subject to 
dynamical Coulomb blockade resemble those of Kondo impurities in Luttinger liquids. Using pre- 
vious knowledge based on the bosonization of quantum impurity models, we show that low- voltage 
conductance anomalies appear at zero temperature. The conductance can vanish at low tempera- 
tures even in the presence of a screened impurity spin. Moreover, the quantitative determination of 
the corresponding Kondo temperature depends on the full frequency-dependent impedance of the 
circuit. This is demonstrated by a weak-coupling calculation in the Kondo interaction, taking into 
account the full distribution P(E) of excited environmental modes. 
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I. INTRODUCTION 

Recent progress in controlling the electronic proper- 
ties of semiconductor-based nanostructures provides new 
ways of probing the physics of strong correlations, and 
envisions a rich interplay with truly mesoscopic effects. 
One of the most prominent examples is certainly the dis- 
covery of the Kondo effect in quantum dots*^ While in 
metals containing a small amount of magnetic impuri- 
ties the increased magnetic scattering of the electrons at 
low temperature results in an increased resistance, the 
mesoscopic realization displays instead a zero-bias peak 
in the conductance^ This setup opens also a way for the 
study of various mesoscopic phenomena in the presence 
of strong correlations, such as nonequilibrium transport 
regimes^ 5 -! 6 - and finite-size effects in the electrodes^ 

One should emphasize that the unitary conductance, 
a feature in the charge sector, is currently taken as the 
fingerprint of the Kondo effect, a property of the im- 
purity spin being screened by the conduction electron 
spins. Yet, it is not obvious that unitary conductance 
and spin screening should always coincide. In fact, one 
should bear in mind that conductance through a single 
or a double junction is sensitive to the circuit environ- 
ment, made for instance of a resistance R. The parameter 
r = R/Rk defines the normalized resistance with Rk = 
h/e 2 . For such an ohmic bath, in a strongly resistive sin- 
gle junction of resistance R t >> Rk and capacitance C, 
with charging energy E c = e 2 /2C, a low-bias anomaly 
dl/dV ~ (V/E c ) 2r appears. This phenomenon is called 
dynamical Coulomb blockade (DCB), (see Refs. (9llToj 
for a review). It has been extensively studied in normal 
single tunnel junctions with resistive leads both experi- 
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well as for general scatterer o 23 ' 24 , in normal double tun- 
nel junctions ; 25 ' 26 superconducting junction s 27 ' 28 ' 29 and 



single barriers in semiconductors^ On the other hand, 
for transparent junctions with R t < Rk, the same phe- 
nomenon occurs, but only below an exponentially small 
DCB energy scale Ec ~ E c e~ RK ' 2Ri . Only when at 
least one channel becomes fully transparent, dynami- 
cal Coulomb blockade disappears and a linear I(V) is 
recovere d 19 ' 23 ' 24 ' 31 ' 32 . Coming back to the Kondo prob- 
lem, one readily sees the difficulty encountered in trying 
to guess the low-temperature behavior. On one hand, 
the zero-temperature regime being fully transparent, one 
might think that the Kondo effect - and the unitary con- 
ductance - are not affected by an ohmic bath. Yet, at 
any finite temperature, the transmission being less than 
one, DCB should reappear. In other words, it is not clear 
whether the Kondo temperature is lower or higher than 
a DCB energy scale. This problem obviously cannot be 
solved without reconsidering how the couplings leading 
to the Kondo effect (and the unitary-conductance peak) 
renormalize at low energy in presence of environmental 
fluctuations. This is the issue we address in the present 
paper. 

Previous investigations of the Kondo effect with dy- 
namical Coulomb blockade focused on the regime where 
the quantum dot is near the charge degeneracy point, 
with a noisy back gate directly coupled to the quantum 
dot j 33 ' 34 ' 35 ' 36 In this case, a Kondo model for the charge 
degree of freedom can be derived, with a direct coupling 
of the dissipation to this charge variable. On general 
grounds this leads to a competition between the Kondo 
screening of the charge doublet by the electrons and the 
localization effect due to the ohmic environment. 

Here, we concentrate rather on the usual spin Kondo 
effect, with environmental electromagnetic fluctuations 
in the electrodes. In this case, dissipation induces a 
markedly different effect on the Kondo physics. Indeed, 
only the inter-lead spin exchange processes involve charge 
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fluctuations across the device and hence dominantly cou- 
ple to the environment. One then assists to a competition 
between several Kondo-type strong-coupling fixed points 
in absence of localization effects. Different conduction 
regimes may then be approached, depending on the dis- 
sipation strength and microscopic couplings. 

At weak dissipation, one expects to preserve the usual 
Kondo effect, but dephasing of the couplings involving 
charge transfer affects the determination of the Kondo 
temperature Tr- below which the spin is eventually 
screened on a quantitative level. For this purpose we ex- 
amine how the perturbative regime of the Kondo model is 
modified by a general dissipative environment in the for- 
malism of the P{E) theory. Focusing on a circuit with 
a finite zero- frequency impedance Z(u> = 0) = R, but 
keeping a complete description of the environment char- 
acterized by the above low-energy ohmic behavior and 
a higher energy tail related to the circuit capacitances, 
we explicitely calculate how those parameters influence 
the temperature scale Tk at which the Kondo resonance 
forms. At temperatures lower than Tk, we can apply 
previous knowledg o 37 i 38 ' 39 i 40 ' 41 on the Kondo effect be- 
tween Luttinger liquids (to which our setup is equivalent 
in some low-energy limit only) to capture the current- 
voltage relation. This extends the previous mappings for 
single or double junction s 21 ' 31 . For the weak-dissipation 
regime, we thus find that nonohmic behavior, thus DCB, 
prevails in presence of particle-hole asymmetry. Particle- 
hole symmetry can be restored for some specific values of 
the dot gate voltage allowing for a unitary conductance, 
albeit with anomalous corrections. 

At large dissipation, which may be more challeng- 
ing to realize experimentally, the electron transfer pro- 
cesses are even more strongly suppressed. As soon as the 
environmental impedance R reaches half the quantum 
value h/e 2 , a generic nonohmic transport regime develops 
(valid independently of particle-hole symmetry, but with 
unbalanced left /right couplings), where the Kondo ef- 
fect occurs through the strongest coupled electrode only. 
With balanced couplings, a two-channel Kondo regime 
may be feasible, and represents a generalization of a re- 
cently proposed 42 ' 43 and experimentally studied setup 4 ^ 
with strong Coulomb blockaded leads, formally realizing 
the limit r ^> 1. We note that our proposal does not 
suffer from a low-energy cutoff set by the level spacing in 
the leads. Notice that weak dissipation (R <C Rk) is not 
sufficient to reach such non- Fermi liquid fixed point. 

The paper is organized as follows. In Sec. [IT] we give a 
characterization of the ohmic environment, and introduce 
the various physically relevant scales for the problem. In 
Sec. IIIII we propose a phenomenological model describ- 
ing how the Kondo effect may be affected by those envi- 
ronmental fluctuations. This model is analyzed both at 
weak coupling, taking into account the full spectral func- 
tion P(E) of the environment, and at strong coupling 
(zero temperature), using a low-energy mapping onto a 
Luttinger chain through bosonization, explicitly derived 
in Appendix |A"1 We conclude the paper with an outlook 



on possible future theoretical developments and experi- 
mental signatures of dynamical Coulomb blockade on the 
Kondo physics. 



II. ENVIRONMENTAL COULOMB BLOCKADE 

A. Circuit theory 

In this part we briefly summarize the basics of circuit 
theory? 9 -*^ This enables us to introduce also some nota- 
tions that will be used along the paper. The device we are 
interested in consists of a nano-object, which can be an 
artificial atom, a molecule or another interacting system, 
connected to reservoirs by two tunnel junctions. Let us 
call S this interacting system. Each time an electron tun- 
nels in or out the system S, this electron can excite modes 
in the electromagnetic environment, describing the cir- 
cuit external to the tunnel junction. As the environ- 
ment comprises S itself, the crucial part is to properly 
describe the tunnel junctions. A semi-phenomenological 
approach able to capture these effects has been developed 
by Nazaro v 10 ! 14 and Devoret et al.M- The environment 
is modeled by its own impedance Z(ui) (now arbitrary), 
leading to an impedance Z t (uS) seen by the junction. We 
denote by Q the charge displacement at the surface of 
the capacitance caused by a tunneling event. Note that 
Q is a collective variable. We also define ip the phase 
conjugate to Q satisfying [<p, Q] = ie. Let us now derive 
the correlation functions of the phase variable ip in the 
case of ohmic fluctuations. 



B. Ohmic phase fluctuation spectrum 

Working with imaginary frequencies, we can 
parametrize the bosonic Green function 45 for the 
phase variable as 



G v (iuj) 



(ip(iu)ip(-iu)) 
Rk i i , 



2n 



i -1 



Z(iu>) 



RkCuj 2 



(1) 
(2) 



where the first contribution in Eq. ^ defines the the 
frequency-dependent impedance Z(uj) of the circuit, 10 
while the second is associated with the junction capac- 
itance C (see Sec. IIII Al for details). We take H = 1 
in the following. For an ohmic environment, the cir- 
cuit resistance is given by R = Z(lj = 0), which from 
Eq. ([2]) yields the natural cutoff frequency lor — 1/(RC). 
However, even for strong dissipation (R ~ Rk) &r is a 
very large frequency, at which Z(y>) cannot be reduced 
to its zero-frequency value. The circuit impedance typ- 
ically decays as Z(u) ~ y/R/\cj\Ci = R^uji/uj for uj 
larger than u>f, u>i = \/(RCi), and Ci is a measure of the 
leads capacitance. 9 To simplify further calculations, we 
will model 46 this tail of the circuit impedance Z(ui) by 
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a pure capacitive decay Z(cu) ~ Rwifto. This allows to 
rewrite the phase correlations as 



R 

G v {iu) = 2n— — 

Hk 



U) 



(3) 



where uj c — wiojr/(uji + ujr) is an effective cutoff. We 
introduce also the related effective capacitance C c g con- 
trolling the high-frequency tail of the phase fluctuations 
lu c = l/(i?C c ff). Using simple complex analysis, we can 
extract the real-frequency spectral function 



R 

Tm G v (u) = 27T— — 
Rk 



n -1 



which obeys 



G v {iu) = - 



de Jm[G v (e)] 



7T iu) — e 

Fourier transforming in imaginary time r yields 
de 1 



GJr) 



-0*-l 



2m[G v (e)]e 



(4) 



(5) 



(6) 



where f3 = l/(ksT) is the inverse temperature, and per- 
forming the analytic continuation t — > it to real time t, 
we get the usual formula 10 



((<p(t) - v(o)Mo)) 

oo 

did 2TZe[Z t (co)] ( cos(u;t) 



(7) 



R K \tanh(f3u}/2) 



i sin(a;f) 



where Z t (u>) = [iu>C + 1/R] 1 is the impedance seen by 
the junction. Since charge transfer processes are associ- 
ated to the operator e llp (see Sec. lIII A"|) . it is useful to con- 
sider the real-time correlation function (e^*' 1 ^') = 
e^cC*), which by Fourier transformation 

P(E) = — j dte J ^ +lEt 



1 

2^ 



(8) 



describes the probability P(E) that the electron ex- 
changes an energy E with the environment during the 
tunneling process. 

At zero temperature and for long time 



Jtp{t) ~ — 2rln(icj c t) 



(9) 



where r = R/Rk is the dimensionless resistance of the 
circuit, which leads to the low-energy behavior— 



P(E) ~ A 



E 2 



,2r 



for E«u c 



(10) 



with A r 



-2rj 



/r(2r). This power law is however cut 



off for energies larger than w c ; in this range the proba- 
bility decays as 



P(E) ~ 2r 



for E > U! c 



(11) 



We finally mention some useful properties^ of the P(E) 
function: i) detailed balance P(E) — for E < at 
T = 0; ii) normalization J Q °° dE P(E) — 1; iii) integral 
equation 



EP(E) = 2r [ duj - 
Jo i 



P(u>) 



j-E 



(12) 



In the next section we analyze how the Kondo effect 
occurring in mesoscopic quantum dots is modified by the 
presence of such an ohmic environment. 



III. THE NOISY KONDO EFFECT 
A. Derivation of the model 

In this section, we consider a single quantum dot 
connected to two metallic leads in the Kondo regime, 
and investigate how the Kondo effect is affected by the 
presence of environmental Coulomb blockade. We as- 
sume that electromagnetic fluctuations are modeled by 
an impedance Z(uj) ~ R closing the circuit as described 
above. The equivalent circuit depicted in Fig. Q] is the 
one of a single-electron transistor. 



dot 



® 



FIG. 1: (Color online) Schematic representation of a circuit 
with a dot coupled to two leads and an impedance Z. 

A model Hamiltonian describing this situation can be 
written as H — H dot + -ffi cads + H tun + ffbath, where 



Hdot = ^ edd v d °- + u nd > 



„t 



(13) 
(14) 



leads — e k C k,a,a C k,a,a J 

fc,(T,a— L I R 

fftun = ^(iLe^dtc^+^e^c^rfJ+H.c. 

(15 



-Hbath describes the environmental degrees of freedom 
consisting of a collection of harmonic oscillators, alter- 
natively expressed by a phase correlator such as Eq. Q . 
In these equations, c k creates an electron with energy 

ek in lead a — L, R with spin a, and d\ creates an elec- 
tron in the dot with spin a. H^ot is the usual Anderson 
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Hamiltonian, where e<j is the dot energy controlled by a 
gate voltage V g , while U describes the Coulomb inter- 
action within the dot. The phases tp a appearing in the 
tunneling Hamiltonian are related to the voltage fluctu- 
ations SV a felt by an electron during a tunneling event 
through the L/R junction as follows 



phase-dependent couplings 



cp a (t) =e SV a (t')dt' 



(16) 



The phases ip a are conjugate to the charge Q a on the 
junction capacitance a such that [ip a , Qa] — is. As they 
originate from the same bath, the phases tp a are clearly 
not independent. In order to relate the voltage fluctu- 
ations in the junction a to Z(u>), we use the standard 
effective-circuit description^ when an electron tunnels 
in the left junction, it feels an environment that can be 
modeled by the right junction in parallel with Z(lo). The 
left/right junctions are simply modeled by a tunneling re- 
sistance R t in parallel with a capacitance Cl/Cr. This 
description is consistent with the single-electron tunnel- 
ing terms in Eq. (|15[) . and may be derived more micro- 
scopically by coupling independently each of the tunnel- 
ing operators to the bosonic charge modes in the cir- 
cuit—^. Neglecting the Kondo correlations for a mo- 
ment, the setup of Fig. [1] is simply a double junction, 
which has been treated in great detail in Ref. [lCj in the 
sequential regime. For a quantum mechanical descrip- 
tion of the double junction in series, it turns out to be 
more convenient to work with the charges q — Ql — Qr 
and Q = C "c^+c^ R being uncoupled in the absence 
of tunneling; q corresponds to the dot charge quantized 
in units of e and Q is simply the total charge carried by 
the capacitance C — q^+q r that couples directly to the 
impedance Z . One may then introduce the phases ip and 
ip conjugate to Q and q respectively, related to the phases 
H> a by ipl = {C/C L )(p + ij} and tp H = -(C/C R )<p + ip. In 
the equivalent circuit, we neglected the part correspond- 
ing to the gate voltage controlling the number of elec- 
trons on the dot. This approximation is justified, as the 
capacitance C g associated to the gate in general satis- 
fies C g <C C a and gate voltage fluctuations usually play 
a minor role (their contribution has been addressed in 
Ref. Si). 

Since we are interested in the Kondo regime, where 
the real charge fluctuations in the dot are suppressed, it 
is convenient to perform a generalized Schrieffer- Wolff 
transformation, or equivalently integrate out the dot 
charge degrees of freedom. The complication due to the 
time-dependence of the tunneling term in Eq. (jT5]) can 
be circumvented by a few additional assumptions. First 
we will assume that \SV\ <C \£d\,U + Ed which consti- 
tutes no serious limitation. We further suppose that the 
energy excitations of the environmental bosonic modes 
lo are small enough to satisfy u <C \sd\,U + Ed- Un- 
der these two conditions, the dot is modeled by a spin 
S described by a Kondo- like Hamiltonian Hk, involving 



Hk = ^ J a ,p < 



a, 



where 



Ef T a,a' <=? 
c k,<j,a 2 k',<r',0 ' ' 

(17) 



k,a,k f a 



Ja.j3 ~ 2t Q t/3 



-1 1 



Id U + £ d 



(18) 



In Eq. (|17p we should a priori include a dressed potential 
scattering term 



a.B 



a,P e L fe,o-,Q L fc' ,a' ,/3' 

k,(7,k' a' 



(19) 



where \ \ 



1t n ti 



(=1 



Nevertheless, the cou- 



LaL P\e d U+e dj 
pling constants V a p turn out to be either irrelevant or 
marginal. As they are not modifying our weak-coupling 
analysis, we neglect these terms. However, these terms 
associated to particle-hole symmetry breaking are impor- 
tant near the strong-coupling fixed point as we will see. 

The omission of the environment fluctuations with re- 
spect to €d or U in Eq. (j 1 8|) means that during a vir- 
tual charge fluctuation or cotunneling event^ no energy 
is exchanged with the environment. This assumption is 
justified by the low value of the circuit cutoff frequency 
Lj[. It contrasts for instance with previous treatments of 
so-called inelastic cotunneling^, or coupling with local 
phonon modest. We expect that going beyond this ap- 
proximation would lead to small corrections to the model 
defined above, and their treatment (see Ref. [5(| for a dif- 
ferent setup and simpler environment) is clearly beyond 
the scope of the present paper. Therefore, the elemen- 
tary Kondo couplings result from quasi-" elastic" cotun- 
neling (in the sense of Ref. 49]), where the initial and 
final states of the spin-flip transition may however in- 
volve different environmental energies. The inter-lead or 
" backscattering" Kondo couplings Jlr and Jrl now in- 
clude the effects of the environment, embodied by the 
dynamical phase <^i(t) — fR{t) = ip(t). This dynamical 
phase clearly modifies the behavior of these couplings, 
which we will analyze in the next two subsections. On 
the contrary, within our quasi-elastic approximation, the 
intra-lead couplings Jll and Jrr are not dressed by 
phase fluctuations. 



B. Weak Kondo coupling analysis 

1. Poor-man's scaling 

The Hamiltonian (fT?]) can be regarded as the usual 
Kondo model, where the inter-lead spin processes J a ,/3, 
with a ^ P transferring charge between the electrodes, 
are dressed by fluctuating phases described by the ohmic 
spectrum of Eq. © . The case of strong Coulomb block- 
ade, corresponding to R = oo, has been previously stud- 
ied in Ref. [43| . The calculation of observables to lowest 
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order amounts to compute the diagrammatics at order 
J 2 . Clearly the renormalization of the J^ R vertex de- 
pends on the combination Jlr(Jll + Jrr), such that a 
single unpaired e lip term occurs, and does not affect the 
usual renormalization equation 



f A 

J'lr. — Jlr + Jlr{Jll + Jrr) / de 

J ASA 



Po(e) 



(20) 



where po(e) — 1/2W is the free electron density of states 
and W is the electronic bandwidth. Following the poor- 
man's scaling philosophy, we introduce here a small shift 
5A in the running cutoff A and computed the result- 
ing change in the Kondo coupling. This flow starts at 
some high-energy initial cutoff D ~ min(U,W), and re- 
duces towards lower energies as long as it is perturba- 
tively controlled. Introducing dimensionlcss couplings 
3a, = PoJ a ,p yields 



dj 



LR 



(flog A 



-JlrUll + Jrr) 



(21) 



On the contrary, the Jll process is obtained from the 
combination J 2 L + J\ R , so that the second term involves 
the combination of two exponential phase factors. This 
amounts to replace in the diagrammatics the free elec- 
tronic Green function Gq(t) = (c r (t)c r (0)) by a mixed 
one: 

G mlx (r) = <4(T)c fl (0)e^) e -^°>) (22) 
« G (r)e J ^ . (23) 

We denote p m ix its associated density of states. Per- 
forming the analytic continuation of the above Green 
function at zero temperature, we find the simple formula 



/•M 

Pmix(u) = P0 dE P(E) 

Jo 



The renormalization equation thus reads 



(24) 



J'll — Jll + Jll 



de PM + jl R r de Pmix(£) 



A-SA ~ e 

(25) 



'A-5A ~ e 
leading to the flow equation 

^- Al - Ar J\eP (E) . (26, 
In the following we analyze this result. 



2. General considerations 



The flow equations (j21j) and (|26[) are the main result of 
this Section. They describe the weak-coupling behavior 
of the Kondo model in a generic environment charac- 
terized by the P{E) distribution, under the quasi-elastic 



approximation. We now specialize to the case of ohmic 
dissipation, for which P{E) was discussed in Sec. HTl 

Although the complete form of P(E) is complicated, 
under the assumption that to c <C U (consistent with the 
quasi-elastic approximation) we can roughly sketch the 
behavior of the flow in a simple two-step picture with re- 
spect to the energy scale ui c . At the beginning of the flow, 
A 3> uj c and the energy integral in Eq. (f2"6"| is unity since 
P(E) is a probability. One thus recovers the usual Kondo 
flow equations, i.e. dissipation is ineffective in this case. 
This can be understood by the disappearance of DCB at 
high energy, and will give rise to an initial strong renor- 
malization of all Kondo couplings as if Coulomb blockade 
was absent. 

This flow continues until A reaches u) c , where the 
anomalous low-energy behavior olP{E) given by Eq. (I10p 
comes into play. At this stage, it is useful t o introduce 
new dimensionless couplings X a ,p = \J A r /2r(A/uj c ) r j a ^ 
for a ^ /3, and X a . a = j a ,a, which allow to rewrite the 
flow equations as 



dX 



LR 



d\og A 
dX LL 
dlog A 



rX LR - X LR {X LL + Xrr) (27a) 

(27b) 



X 2 X 2 
~ A LL — A LR 



Let us first give some general comments on these equa- 
tions. First we notice that the dissipation affects the 
inter-lead couplings Eq. (|27a|) and drives them irrele- 
vant, at least as long as Xll + Xrr is not too large. 
This is consistent with the picture where charge transfer 
across the quantum dot is suppressed by DCB, and has 
the clear consequence of diminishing the Kondo temper- 
ature at increasing values of r (see also the discussion 
in Sec. IIII B~3l and corresponding Fig. [J). Second, these 
flow equations are reminiscent of those obtained for the 
tunneling through a magnetic impurity in a Luttinger 
liquid, 38 although we did not obtain them from bosoniza- 
tion. This identification will be put on rigorous grounds 
in Sec. IIII CI by virtue of an exact mapping between the 
two problems. This also allows the computation of trans- 
port properties of the device at low temperature, that the 
weak-coupling calculation is unable to capture. 

The above two-step argument thus shows that the slow 
onset of DCB at the scale oj c reinforces the flow with 
respect to a pure Luttinger flow of the type (|27p. so that 
the decrease of the Kondo temperature with increasing 
dissipation will not be as pronounced as in the case of a 
pure Luttinger chain. 

We also note that the cutting of P(E) at ui c is certainly 
not quantitative since P(E) is a broad distribution, and 
we will solve the flow Eqs. (|2"Tj) and (|2"6")l numerivally, 
with the full evaluation of P(E) according to the integral 
Eq. lfl"2]) , as performed below. 

A final remark concerns the case of strong Coulomb 
blockade for r>l. In this regime, the Coulomb charging 
energy E c = e 2 /2C emerges as a natural cutoff 56 with 
P(E) ~ 6(E — E c ). For this case the two-step argument 
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becomes exact^ and leads to a strong renormalization 
of the Kondo temperature. 



dependent impedance of the junction in the circuit) af- 
fects the determination of the Kondo temperature in a 
significant quantitative manner. This is also seen by a 



3. Kondo temperature 



We start for simplicity by considering the solution of 
the pure Luttinger liquid case, described by Eqs. (f!?T)) . 
This corresponds to the regime of large w c , ignoring the 
mixing of boson excitation energies in the initial Kondo 
coupling. The results displayed in Fig. [5] however illus- 
trate in the simplest manner the role of the parameter 
r. The first general effect is, as expected, a systematic 
decrease of the Kondo temperature Tk with increasing r, 
marking a systematic onset of a strong-coupling regime. 
We note that Tk remains nonzero even for strong dissipa- 
tion, contrarily to the noisy "charge" Kondo effect ) 33 i 34 in 
which voltage fluctuations act as a fluctuating magnetic 
field would act in the spin Kondo problem, that may 
prevent the formation of the Kondo resonance. Apart 
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FIG. 3: (Color online) Flow of the Kondo coupling jt,jj(A/D) 
(full lines) according to Eqs. (|2ip and (|26[1 . and the corre- 
sponding integrated distribution J A dE P(E) (dashed lines), 
for the dissipation strength r — 1, initial couplings j'll = 
Jlr =0.1 and several values of the dimensionless effective 
capacitance c c s — DRkC c s = 1,2,5,10,20 (from bottom to 
top). 



systematic study of the Kondo temperature (defined here 
by the criterion Jll(^- — Tk) — 10) as a function of dissi- 
pation r and for several values of the effective impedance 
C cff in Fig. m 



FIG. 2: (Color online) Flow of the Kondo couplings 
Xlr(A/ D) (full line) and All (A/7?) (dashed line) according 
to Eqs. (|27p (i.e. pure Luttinger liquids), for different values 
of r. The initial values of the Kondo couplings at scale A = D 
are taken here as All = Alb = 0.1 (with L/R symmetry). 

from this quantitative effect, we infer from Fig. fj] that 
dissipation can furthermore discriminate between several 
strong-coupling fixed points. Indeed when the impedance 
R becomes comparable to the quantum value Rk, the 
inter- lead couplings Xlr are strongly driven to zero. This 
should clearly lead to qualitatively different transport 
properties compared to the case of weak dissipation. This 
also allows for a stabilization of a two-channel Kondo 
effect by DCB, for the case of balanced All and Ark 
couplings. These aspects are confirmed by the strong- 
coupling analysis performed in Sec. IIII CI 

The solution of the flow Eqs. (p?Tj) and ptj|) with ohmic 
dissipation, using the full P(E) function determined by 
Eq. (H2), is shown in Fig.[3J It illustrates the general fact 
that the full distribution function P(E) of the environ- 
mental modes (or equivalently the complete frequency- 




0.2 0.4 0.6 0.8 



r 

FIG. 4: (Color online) Kondo temperature Tk as a function of 
dissipation r, for several values of the dimensionless effective 
capacitance c c g = 1, 2, 5, 10, 20 (from bottom to top). 



7 



C. Transport anomalies 

1. Equivalence to the Kondo model between Luttinger liquid 

leads 

The previous weak-coupling analysis already informs 
about the type of ground state favored by the environ- 
ment. More sophisticated tools are however necessary 
to make precise statements on the nature of the ground 
state, and to compute low-temperature and low-voltage 
properties. Here we will formulate the low-energy equiv- 
alence between the Hamiltonian (jTTJ) and the problem 
of a S = 1/2 magnetic impurity weakly coupled to two 
Luttinger liquids, in the tunneling geometry. This equiv- 
alence has been already shown at the Hamiltonian level 
for a nonmagnetic impurity by Safi and Saleur— (see also 
Ref. [H]]), and can be easily extended to the Kondo model 
through bosonization of the lead electrons^ as derived 
in App. El The final result of this calculation allows to 
identify the effective interaction parameter of the Lut- 
tinger chain as 



K 



1 



l + 2r 



(28) 



The problem of a magnetic impurity in a Luttinger liquid 
has been analyzed using weak-coupling renormalization- 
group equations and a stability analysis of the strong- 
coupling limit by Fabrizio and Gogolin^, see also 
Ref. [39|] for a conformal field theory analysis. It is thus 
not surprising that the low-energy limit (j2"T)) of the flow 
Eqs. (|2"Tj) and (f2T)]) reproduces those results. We empha- 
size that the flow equations obtained for the alternative 
"side" geometr y 51 ' 52 differ from the present "tunneling" 
geometry ones by a relevant scale dimension 1 > K of 
the inter-lead couplings rather than the irrelevant scale 
dimension 1/K > 1 (as r > 0) found above. The result- 
ing physics is thus markedly different for the two cases. 



2. Strong- coupling analysis 

The low-energy physics and therefore transport and 
thermodynamic properties are dominated by the vicinity 
of the strong-coupling fixed point. In order to analyze 
the low-energy physics, two distinctive regimes arise de- 
pending on the value of r. 

Small values r < 1/2. The previous weak-coupling 
analysis has shown that all Kondo interaction terms are 
driven to strong coupling such that the usual one-channel 
Kondo physics should a priori be recovered. However, 
for a more precise statement the vicinity of this fixed 
point should be considered. A detailed analysis shows 
that particle-hole symmetry-breaking terms are relevant 
at the one-channel Kondo fixed point. One indeed finds 
that potential scattering terms like Vc^c in Eq. (TT9)) have 
a scaling dimension (1 + K)/2 = (l + r)/(l + 2r) < 1 and 
are therefore relevant. In a one-dimensional infinite Lut- 
tinger liquid chain, such a term "cuts" the chain into two 



semi-infinite chains below an energy scale T*M. There- 
fore, the ultimate fixed point which the system reaches 
for T = eV = corresponds to an insulating one i 38 ' 39 ' 40 
The analysis of the operator content at this insulating 
fixed point allows to predict the scaling behavior with 
temperature or bias voltage. Here we focus on the dif- 
ferential conductance dl/dV between the left and right 
lead. The leading irrelevant operator corresponds to a 
hopping term between the two leads and has dimension 
(1 + l/K)/2 = 1 + r > 1. We note that particle-hole 
symmetry-breaking effects can be tuned to zero, using 
the dot gate voltage. The Kondo resonance may therefore 
develop in transport for some particular value of the gate 
voltage V* , as specified in Eq. (|29b[) below. The finite- 
temperature corrections to the unitary limit are given 
by the leading irrelevant operator nearby this perfectly 
transmitting fixed point when particle-hole symmetry is 
imposed. This operator has dimension 2K and is irrel- 
evant for K > 1/2i2S These results imply an anomalous 
behavior of the differential conductance at low tempera- 
ture: 



dl_ 

dV 

dl_ 

dV 



(29a) 



— — ~ 2G 



K 



2r 



iovV g = V* (29b) 



with Gk = 1/Rk — e 2 /h, and a, b dimensionless nonuni- 
versal prefactors. At finite temperature T and zero 
bias voltage, the energy eV should be replaced by T in 
Eqs. (|29[). The conductance in Eq. (|29a[) is similar to the 
conductance through a single tunnel junction. One may 
also wonder how the conductance is modified when V g 
slightly deviates from V*. Since particle- hole symmetry 
is relevant at this perfectly transmitting fixed point, and 
the resonance is very sharp. Following Ref . [37j , it can be 
shown that the width of the resonance scales as T^-~ K " 2 
where V g deviates from V* . The differential conductance 
escapes from the unitary limit as follows 



dl 

dV 



2G K 



eV 



for v g « v; 



(30) 



where T is a characteristic energy scale related to SV g 



V„ — V* , which vanishes as SV, 



2/(l-K) 



at resonance^ 



'g ' g ) .nmoii« 0,0 Wg 

Large values r > 1/2. For large values of r, dissi- 
pation dominates the renormalization-group flow in Eq. 
(|27a[) and the Kondo coupling Jlr renormalizes to zero. 
Therefore coherent charge transfer between left and right 
lead is suppressed, without however completely inhibit- 
ing the Kondo effect, similarly to what occurs in the 
case of strongly Coulomb blockaded lead a 42 ' 43 . Indeed 
the Kondo couplings Jll and Jrr are still driven to 
strong coupling according to Eq. (|27b[) . A two-channel 
Kondo fixed point can be reached provided Jll ~ Jrr- 
By analyzing the operator content of the strong-coupling 
two-channel fixed point, it has been shown to be sta- 
ble against particle-hole symmetry-breaking terms for 
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K < 1/2, or equivalently r > 1/2^.^ Indeed, the 
leading irrelevant operator near the two-channel Kondo 
fixed point occurs in the flavor sector and has dimen- 
sion 1/(2K) = 2r + 1, and is therefore irrelevant for 
r > 1/2^ This operator will therefore dominate the 
transport properties below the two-channel Kondo tem- 
perature Tick- For Jll ^ Jrr the two-channel Kondo 
fixed point is unstable and a single-channel Kondo ef- 
fect occurs in the strongest coupled channel (either left 
or right). The transport for this situation is character- 
ized by the leading irrelevant operator corresponding to a 
hopping term between the two leads as for r < 1/2 with 
dimension (1 + l/K)/2 = 1 + r > 1. From this above 
analysis we can infer the anomalous behavior of the dif- 
ferential conductance dl/dV between the left and right 
leads at low energy: 



dl 

dV 

dl 
dV 



\Twk) 



for Jll ^ Jrr, 



(31a) 



eV 



2CK 



2r-l 



for J LL = J RR , (31b) 



where c, d are nonuniversal numbers and Tick is the one- 
channel (eithr left or right) Kondo temperature. Notice 
the absence of Gk, peculiar to the two-channel Kondo 
effect in Eq. (|31b[) . due to the tunneling geometry. A 
third, nonnoisy electrode would be necessary to obtain 
this behavior, see Ref. [42|j43j]. 

This analysis can be directly extended to a quantum 
dot operating as a spin-1 impurity. The low-energy 
physics resembles then that of a spin-1 impurity embed- 
ded between two Luttinger liquids, recently studied in 
Ref. [H. 



IV. CONCLUSION 

In this paper we investigated the impact of environ- 
mental Coulomb blockade on the Kondo effect in the case 
of moderate to highly resistive leads. An ohmic resistance 
of the environment exceeding i?_R-/2 induces a suppres- 
sion of the inter-lead Kondo interactions, without how- 
ever preventing the formation of a strong-coupling state, 
due to the remaining intra-lead processes. When the tun- 
neling amplitudes between the left and right leads are 
equal, even a two-channel Kondo effect can be reached. 
For an environmental resistance smaller than Rk/2 the 
Kondo effect can fully develop between both leads. How- 
ever, the fully transparent fixed point is stable only 
when particle-hole symmetry is maintained with the dot 
plunger gate voltage. On a qualitative level these results 
imply anomalous low-temperature transport properties 
through the device, even though the Kondo effect sur- 
vives. More quantitatively, the Kondo temperature be- 
low which those features start to appear was shown to 
depend sensitively not only on the dissipation strength, 
but rather on the full spectral function of the excited 
environmental modes. Although strong-resistive effects 



may be difficult to detect in actual experiments, due to 
the difficulty in the realization of well-coupled highly dis- 
sipative setups in semiconducting quantum dots, such a 
reduction of the Kondo temperature at moderate dissi- 
pation as well as nonlincarities in the I(V) characteristics 
should be observed. 

Theoretically, it would be interesting to go beyond the 
quasi-elastic approximation by which our phcnomcnolog- 
ical Kondo model is applicable. This approximation may 
rather underestimate the environmental effects. A gener- 
alization of this study to orbital Kondo effects or mixed 
spin-orbital Kondo setups can also be performed along 
the same line^ merging the present situation with the 
previously studied decoherence effects on the Kondo ef- 
fect in the charge sector. 
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APPENDIX A: EQUIVALENCE BETWEEN THE 
NOISY KONDO PROBLEM AND THE KONDO 
MODEL IN LUTTINGER LIQUIDS 

We provide here the derivation for the low-energy map- 
ping of the Kondo model (|17[) in an ohmic environment 
onto the usual Kondo model between Luttinger liquids. 

The free electronic part of Eq. (TT1)) is first conveniently 
expressed as a chiral Hamiltonian for left and right elec- 
trodes 



Pleads — —ivp 



/+oo 
-OO „ n 



(Al) 



We bosonize each channel, keeping the Klein factors ex- 
plicitly 



^=b aa exp i — 



V2 



(A2) 



This allows to express the above Eq. (|Alj) in terms of 
charge and spin bosons in each channel 

r+oc 

Pleads = -f dx J2( d ^c a ) 2 + {d x 4>s a f ■ (A3) 



Using the identity for the electron density operator 
1 d x cf> ca + {-l) a d x (t> sa 



c aa C aa 



2tt 



V2 



(A4) 
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and the bosonization dictionary (|A2|) . we can now re- 
express the Kondo Hamiltonian (jTTJ). To clarify the re- 
sulting expression, we define the following basis change 
between the L, R channels 



D C R 



r> sf 



V2 




4>cL - 1 




V2 




4>sL + <? 


hR 


v/2 




4>sL - <? 


hR 



V2 



(A5a) 
(A5b) 
(A5c) 
(A5d) 



We have also to perform a similar transformation on the 
Klein factors^ and introduce for this purpose four alter- 
nate Klein fermions, F c , Ff, F s , F s f, which obey 



We see that the environmental boson ip occurs nat- 
urally as a systematic shift of the "flavor" field cj)f de- 
scribing the charge transfer between the two electrodes. 
To define a simple effective environment, we consider 
the combined effective action for both the local density 
bosons (integrating all x ^ in Eq. (IA3|) ). and the envi- 
ronmental boson from Eq. @ 



So = ^£Kl(Vl 2 + l ( M 



°sf\ 



Rk 
~2R 



M 2 + I^'l 2 ) 



F l f F t 
F^F 



/?f p 
pt p 
/?t p 

r e] r 2] 



(A6a) 
(A6b) 
(A6c) 
(A6d) 



All other bilinears of the original Klein factors are eas- 
ily computed from the above expression, e.g. F^F 21 - — 



pf p pt p 



p p pt pt 



F.Fl 



We introduce also spin anisotropies of the exchange 
constants, with the new notation (Jp, Jp, Jg, J|), 
where the subscript F (B) stands for Jll = Jrjj (Jlr = 
Jrl), and the superscript _L (z) for the x, y (z) spin orien- 
tations respectively. We assumed here perfect symmetry 
between the left and right couplings not to overburden 
the notation, although this will not affect the validity of 
the mapping we want to establish. The final expression 
for the Kondo term (flT)) then reads 



where 



V = 




2R, R K 1 



(A8) 
(A9) 



and 



K 



1 



1 + 2R/R K 



(A10) 



H 



K 



pl 

Aira 

—S z d x <f) 
4ira 



i<j> B f 



B S'Fte 



i(j>f —itp 



t p -i<t>sf 



H.c. 



Fje 



i(f> f —iip 



-i(fi s 



P P * 



H.c. 

f H.c. . 

(A7) 



Therefore the role of the ohmic environment on the 
Kondo effect is to replace all exponential terms e^f in 
the bosonized expression by ^tl^^ . Starting from a 
standard magnetic impurity between two Luttinger liq- 
uid leads characterized by an interaction parameter K 
(without environmental modes) , with the open boundary 
bosonization condition^ leads to the same final expres- 
sion, thus proving the desired equivalence. 
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